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An improved model for the transit entropy of monatomic liquids 
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In the original formulation of vibration-transit (V-T) theory for monatomic liquid dynamics, the 
transit contribution to entropy was taken to be a universal constant, calibrated to the constant- 
volume entropy of melting. This model suffers two deficiencies: (a) it does not account for exper- 
imental entropy differences of ±2% among elemental liquids, and (b) it implies a value of zero for 
the transit contribution to internal energy. The purpose of this paper is to correct these deficiencies. 
To this end, the V-T equation for entropy is fitted to an overall accuracy of ±0.1% to the available 
experimental high temperature entropy data for elemental liquids. The theory contains two nuclear 
motion contributions: (a) the dominant vibrational contribution S V ib{T/9o), where T is temperature 
and 6*o is the vibrational characteristic temperature, and (b) the transit contribution Str{T/9 tr ), 
where 0t r is a scaling temperature for each liquid. The appearance of a common functional form of 
Str for all the liquids studied is a property of the experimental data, when analyzed via the V-T for- 
mula. The resulting Str implies the correct transit contribution to internal energy. The theoretical 
entropy of melting is derived, in a single formula applying to normal and anomalous melting alike. 
An ab initio calculation of 8o, based on density functional theory, is reported for liquid Na and Cu. 
Comparison of these calculations with the above analysis of experimental entropy data provides 
verification of V-T theory. In view of the present results, techniques currently being applied in 
ab initio simulations of liquid properties can be employed to advantage in the further testing and 
development of V-T theory. 

PACS numbers: 05.70.Ce, 61.20.Gy, 64.70.dm, 71. 15. Mb 



I. INTRODUCTION 

We are interested in the theoretical description of the 
motion of nuclei (or atoms) in real monatomic liquids. 
For many years, such descriptions have been available for 
gases and crystals, consisting in each case of an approx- 
imate but tractable "zeroth order" Hamiltonian, plus 
complicated but small corrections. Zeroth order for a gas 
is free particle motion, from Boltzmann [1| , and the cor- 
rection is potential energy. Zeroth order for a crystal is 
harmonic vibrational motion and the major correc- 

tion is anharmonicity. These theories are extremely valu- 
able, as the zeroth order Hamiltonian provides a complete 
orthogonal basis set for the nuclear motion. Hence for 
any physically meaningful problem, the motion can be 
analyzed and statistical mechanics can be constructed in 
terms of the basis set. These theories account for equilib- 
rium and noncquilibrium properties of gases and crystals 
to an accuracy on the order of experimental accuracy. 

Our proposal for this kind of theory for monatomic 
liquids is vibration-transit (V-T) theory The key 
postulate is that the many-body potential energy surface 
is overwhelmingly dominated by intersecting macroscopi- 
cally equivalent random valleys. The zeroth order Hamil- 
tonian expresses normal mode vibrational motion in a 
single random valley harmonically extended to infinity. 
"Macroscopic equivalence" means that this Hamiltonian 
for any random valley gives the same statistical averages 
in the thermodynamic limit [a]. The motion of nuclei is 
then composed of two parts: brief periods of vibration in 
one random valley, interspersed with transits which carry 
the system between neighboring random valleys. Tran- 
sits, each involving a small local group of nuclei, proceed 



at a high rate throughout the liquid. For calculation of 
the partition function, the effect of transits is to correct 
the potential surface for the valley-valley intersections 
For nonequilibrium calculations as, e.g., of time 
correlation functions, the same transits provide the diffu- 
sive jumps of the nuclei 0,0|- The vibrational motion is 
tractable, and is calibrated from potential properties of 
a single random valley. The vibrational contribution to a 
thermodynamic function is around 90% of the total 0, [f| . 
The transit motion is complicated, but its contribution 
to a thermodynamic function is only around 10% [J, [y]. 
This paper is mainly concerned with the vibrational and 
transit contributions to the entropy of monatomic liquids. 

In the original formulation of V-T theory, transits are 
accounted for only insofar as they give the liquid access 
to all the random valleys. This property is modeled by 
multiplying the single random valley partition function 
by a universal number, calibrated from entropy of melt- 
ing data [H . For St r , the transit contribution to entropy, 
this yields St r — OSNks, the same for every monatomic 
liquid. Theory for the total entropy agrees with high- 
temperature experimental data for normal melting ele- 
ments to within ±0.2iVfc B ( 4], Fig. 2). Since 0.2Nk B is 
approximately 2% of the total entropy for the elemental 
liquids, the original V-T formulation is quite satisfactory 
for such a simple model. Nevertheless, the original for- 
mulation suffers two deficiencies which we wish to correct 
here. 

1. Since St r is a universal constant, it does not ac- 
count for the different behaviors of individual liq- 
uids. However, these differences largely account for 
the scatter which results in theoretical errors of up 
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to 2% in the entropy. 



The characteristic temperature 9q is given by 



2. Since S tr does not depend on V or T, the consistent 
contribution to all other thermodynamic functions 
is zero. But we now know that the transit contri- 
bution to internal energy at melt is around 10%, 
and this energy must be included to get a good 
theoretical value of the melting temperature T m . 

To correct these deficiencies, we must account explicitly 
for the V and T dependences of Str- 

To this end, in Sec. |TT] we extend the original analysis 
of experimental high-temperature entropy data to nor- 
mal and anomalous melting elemental liquids. We find 
that all the data for St r (T) at constant volume can be 
fitted to a single curve, providing a scaling formula for 
the V and T dependences of S tr - The fitting yields values 
for each liquid of the vibrational characteristic tempera- 
ture #o and of a new transit characteristic temperature 
0t r . The original "universal entropy constant" has a role 
in the present formulation, where it is expected to vary 
weakly with V, and also to vary weakly among the ele- 
mental liquids. 

In liquids, as in all condensed matter systems, the po- 
tential energy that governs the nuclear motion is the elec- 
tronic ground state energy as a function of nuclear posi- 
tions 1 2, 6]. Hence the parameters in the liquid dynamics 
Hamiltonian can be calculated from electronic structure 
theory. Over the years, pseudopotential perturbation 
theory for nearly-free-electron (NFE) liquids has been ex- 
tremely useful in the study of liquid dynamics for real liq- 
uids 

SHIM EH, 

[l2l ] . This is the basis of a series of tests 
of V-T theory for liquid Na @, El El EE El • 

Extension 

of this principle beyond NFE liquids, by means of density 
functional theory (DFT), is the subject of Sec. 11111 First, 
V-T theory for the entropy of melting is reviewed, and a 
single equation covering normal and anomalous melting 
is derived. Results from a new method [It], El], using 
DFT to calculate the vibrational parameters, are then 
applied to test the V-T theory of entropy for Na and 
Cu. Ab initio techniques are currently being applied to 
a wide range of liquid dynamics studies, and their poten- 
tial in testing and developing V-T theory is noted. In 
Sec. IIV1 broader application of the present reformulation 
of the transit entropy is described. The verification of 
V-T theory provided by the present ab initio test is also 
discussed. 



\n(k B o ) = (ln(M), 



(3) 



where (...) is the average over the vibrational normal 
mode frequencies u>. The quantum corrections to Eq. ^ 
are straightforward and are negligible in the present 
analysis. Str represents the transit motion of the nu- 
clei. S e i represents the excitation of electrons from their 
ground state with nuclear positions fixed at a random 
structure. Two small contributions, neglected here, ex- 
press anharmonicity of the vibrational motion and the 
interaction between nuclear motion and electronic exci- 
tations (@, Sec. 4). 

The experimental data are at ambient pressure, where 
the volume increases with temperature. It is most help- 
ful to remove the volume dependence of the experimental 
data, by correcting S expt (V, T) to S expt {V m ,T), where V m 
is the fixed volume of the liquid at melt. At fixed vol- 
ume, the parameters of our analysis are simply constants. 
With (dS/dV)T — (3Bt, where (3 is the thermal expan- 
sion coefficient and Bt is the isothermal bulk modulus, 
the correction to second order is 



S e xpt\Ymt T) 



1 



S expt (V,T)+r)VpB T +^r, 2 V 2 . 



f d((3B T y 
dV /T 
(4) 

where 77 = (V m — V)/V. For a given liquid the anal- 
ysis requires highly accurate experimental data for the 
entropy and its first volume correction, at a significant 
range of temperatures above T m . The liquids satisfying 
this requirement are ten NFE metals. These are listed in 
Table [U along with T m , the highest temperature of our 
analysis Th, and references for the experimental data. 

The second volume correction is negligible until the 
first volume correction reaches a magnitude around 
O.SNks- The second volume correction is calculated for 
Hg, since the data are sufficient for this purpose [Hj]. For 
the alkali metals, the second volume correction is esti- 
mated from experimental data in the vicinity of the melt 
curve 34]. For the remaining five liquids, the highest 
temperature Th is such that the second volume correc- 
tion can be neglected. 

The analysis will now be confined to the volume V m , 
and the corresponding notation will be suppressed. The 
first step is to find St r from S exp t, using Eq. (JTJ) in the 
form 



II. ANALYSIS OF EXPERIMENTAL ENTROPY 
DATA 

In V-T theory, the liquid entropy is given by 

S(V, T) = S mb (V, T) + S tr (V, T) + S el (V, T). (1) 

S v ib describes the nuclear motion in a single random val- 
ley harmonically extended to infinity. In classical statis- 
tical mechanics, 



Svib(V, T) = 3Nk B {In [T/0 o (V)] + 1} . 



(2) 



Str(T) = S expt (T) — S v ib(T) — S e l(T). 



(5) 



SeiiT) is calculated from free electron theory in the 
leading Sommerfeld expansion. This is sufficiently ac- 
curate for the liquids studied here, because S e i(T) < 
0.02S expt {T) throughout the analysis. S V ib(T) is calcu- 
lated from Eq. © with an initial choice for 9q. At this 
point, eight curves of St r vs. T have clear maxima, with 
the exceptions being Pb and Ga. For the eight, the tem- 
perature at the maximum is denoted 9 tri and we graph 
Str(T) vs. T/0tr- The curves look like they will fall on 
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TABLE I: Results of the high-temperature entropy analysis for ten liquids. T m is the melting temperature and Th is the highest 
temperature of the analysis. Experimental data are from the references cited. Sn and Ga are anomalous melters [ljj. S nuc is 
given by Eq. (|6]), and (5S exp t) is the estimated mean error in the high-T entropy data after correction to V = V m (see Appendix 
III). 



Element 


T m (K) 


T h (K) 


References 


ft, (K) 


6tr (K) 


S nU c{T m ) (Nk B ) 


(5S expt ) (Nk B ) 


Na 


371 

Ol l.U 


1 1 00 

1 1UU 


[on 91 99 9"i\ 

L£!i' zl - zz - ZJ| 


Q7 fi 


^70 

O i yJ 


7 79^ 


07 


K 


336.4 


1040 


[22, 23, 24, 25] 


58.0 


570 


8.986 


0.11 


Rb 


312.6 


900 


[22, 24, 26] 


35.8 


530 


10.183 


0.12 


Cs 


301.6 


948 


[22, 24, 26] 


26.0 


540 


11.032 


0.11 


Al 


933.5 


1400 


[24, 27, 28, 29] 


198.0 


980 


8.451 


0.09 


Pb 


600.6 


1023 


[21, 27, 28, 30] 


53.3 


580 


11.041 


0.09 


In 


429.8 


920 


[21, 31J 


74.2 


600 


9.040 


0.16 


Hg 


234.3 


630 


[21, 32] 


52.5 


260 


8.284 


0.04 


Sn 


505.1 


1173 


[21, 27, 30] 


73.7 


640 


9.567 


0.15 


Ga 


302.9 


773 


[21, 33] 


99.6 


360 


7.148 


0.07 



a single line if they are shifted by various constants in 
Str- This is done by varying 6$, since a change in 9q 
changes S tr by a constant [Eqs. ^ and |J5J)]. To bring 
the curves together, a common value for the maximum of 
Str is needed. We choose 0.8Nk B , the universal entropy 
constant of the original V-T formulation. All ten liquids, 
including Pb and Ga, can be shifted to lie on a single 
smooth curve, as shown in Fig. [TJ 



at T — 9 tr - For each liquid, this information plus the 
parameters 0q and 9t r constitute a fit of the experimental 
entropy data at volume V m . However, we do not expect 
a precise common maximum to hold for all monatomic 
liquids, nor do we expect this maximum to be volume 
independent. For a more quantitative analysis of this 
issue, we turn next to the entropy of melting. 
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FIG. 1: (Color online) Transit entropy for ten liquids cal- 
culated from the data referenced in Table |U Within small 
scatter, all lie on a single curve as a function of the scaled 
temperature, where 9 tr is a material parameter. The scatter 
is small compared to the mean experimental error (Table |l}. 

In this analysis, the actual data for St r {T) for each liq- 
uid is a set of points with small scatter. Therefore, plac- 
ing the ten data sets on a single curve is not a precise 
operation. But the scatter is so small that this impre- 
cision is negligible. The fitted values of 0$ and 9 tr are 
listed in Table HI The total nuclear motion entropy is 



S n uc(T) — S v ib(T) + Str(T). 



(6) 



To show its magnitude, S nuc (T m ) is also listed in Table U 
Aside from experimental data, the information used in 
the analysis is Eqs. (fTJ-fiiJ), plus the constraint that all 
the curves must have a common maximum of 0.8Nk B 



III. TESTING V-T THEORY 

A. Theory for the entropy of melting 

The experimental information relating to the univer- 
sality of the transit entropy is the entropy of melting. 
To use this information we need to express the entropy 
of melting theory in terms of the entropy formulation of 
Sec. |TTJ To keep the notation simple we shall continue 
with our standard notation for the liquid, and use a su- 
perscript c to denote the crystal. Hence S and S c are the 
entropy of the liquid and crystal respectively, and V rn is 
the liquid volume at the melting temperature T rn . The 
measured constant-pressure entropy of melting, corrected 
so that both crystal and liquid have the same volume V m , 
is 



AS(V m ,T m ) - S{V m ,T m ) - S c (V m ,T m ). 



(7) 



For all the normal melting elements for which accurate 
experimental data are available to determine this quan- 
tity, the mean and standard deviations are @, [l9l . l3a ] 



AS(V m ,T m ) = (0.80 ± O.W)Nk B . 



(8) 



The distribution is essentially the same when AS is evalu- 
ated at the crystal volume V^, so the volume dependence 
is weak. Among the elements for which accurate experi- 
mental data are available, six do not belong to this dis- 
tribution and are called anomalous. Their values of AS 
in units of Nk B are 1.48 (Sn), 2.37 (Ga), 2.68 (Sb), 2.62 
(Bi), 3.76 (Si), and 3.85 (Ge). The anomalous AS are not 
only larger than the normal value, they are much larger. 
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As a fiducial for elemental metals, the liquid entropy at 
melt is S m w ION kg. In comparison, experimental error 
is w 0.5%, the width of the normal AS distribution is 
very small at w 1%, and the anomalous A5* is very large 
at 10 - 30%. 

To rationalize these experimental results, we shall 
write the V-T equation for AS(V m ,T m ). The crystal 
entropy, with classical harmonic vibrations, is 

S%V, T) = 3Nk B {In [T/9 C (V)} + 1} + S c el (V, T). (9) 

Then with Eqs. (□} and ©, 

AS(V m ,T m ) = 3Nk B ln[e c (V m )/e (V m )] 

+S tr (V m ,T m ) + AS el {V m ,T m ).{lQ) 

In normal melting, there is no significant change in the 
electronic structure, so that the internuclear forces are 
nearly the same in liquid and crystal, and so is the elec- 
tronic density of states. Hence #o(^m) ~ ^o(Kn) an d 
AS el (V m , T m ) « 0, so that from Eq. (TJJJ, AS(V m ,T m ) « 
S tr (V m ,T m ), It then follows from Eq. © that 

S tr (V m ,T m )^(0.80±0.10)Nk B (11) 

for normal melting elements. 

In contrast, anomalous melting is accompanied by a 
change in the electronic structure. Si and Ge melt from 
covalent crystal to metallic liquid ([3a]. Chap. 3), while 
Sb and Bi melt from semimetal crystal to metallic liquid 
(03> P- 81). For S n and Ga, the electronic structure 
change upon melting becomes apparent through com- 
pression. Si, Ge, Sb, and Bi all have a triple point on 
the melt curve at modest compression (3l}, and so do 
Sn and Ga [13, The interpretation is that com- 

pression drives a relative shifting of electronic bands, 
the shift being continuous with compression for the liq- 
uid, but being concentrated at the crystal-crystal tran- 
sition in the solid. Hence the melting is anomalous in 
the vicinity of a triple point. The classic example is Cs 
[39l liH |4]| . where melting is normal at ambient pres- 
sure, but becomes anomalous under compression when 
the triple point is approached (@|, Fig. 26.5). 

Because of the change in electronic structure, the in- 
ternuclear forces and electronic density of states are sig- 
nificantly different between liquid and crystal. Hence in 
addition to the normal St r ~ 0.8Nk B on the right hand 
side of Eq. ([TO)) , the terms in nuclear vibration and elec- 
tronic excitation are both important. These terms have 
positive sum because melting is entropy driven, and the 
term in ln(0g/0o) is usually dominant. 

Let us apply these results to Fig. [T] Equation 
(fTTj) holds for normal melting, and since St r (V m ,T) 
changes little between T m and 8 tr , according to Fig. [TJ 
Str(V m ,0tr) ~ 0.8Nk B for normal melting elements. 
But this is a purely liquid quantity, independent of the 
melting process, so this relation should be valid for 
monatomic liquids in general. Let us denote this com- 
mon maximum x{V), a function of volume. The above 



TABLE II: DFT calculations of 9 for the liquid at V m HH, 
compared with the same quantity determined from experi- 
mental entropy data. 



Liquid 


P (g/cm 3 ) 


O (DFT) (K) 


o (expt) (K) 


Na 


0.935 


98.4 ± 3.0 


97.6 ± 2.3 


Cu 


8.00 


171 ±5 


171.4 ±5.4 



argument suggests 

X (V m ) « 0.8Nk B . (12) 

In Fig. [1] the approximation is taken to be an equality. 
In principle, x{V) i s a material parameter, but as the 
following test shows, we are not yet able to resolve specific 
material dependence. 

B. Testing through ab initio calculations 

The test reported here is the comparison of 9q from ab 
initio calculations with results from the analysis of exper- 
iment for Na and Cu. The technique [l?], EH calculates 
the electronic ground state by DFT (using the VASP code 
[43 |) for a system of N atoms in a cubic cell, with periodic 
boundary conditions on the nuclear positions. The sys- 
tem is quenched to a structure, where the frequencies and 
eigenvectors of the normal vibrational modes are calcu- 
lated, and 9q is evaluated from Eq. ([3]). The eigenvalues 
(mass times squared frequencies) are always positive, ex- 
cept for the three expressing translation, which are zero 
to numerical accuracy. The structures found are numer- 
ically dominated by random ones. Each calculation is 
done at the density of the liquid at melt, and 6* (DFT) is 
listed in Table HU 

Table [TT1 also lists $o(expt), which is from Table Q] for 
Na. For Cu, with insufficient data for the analysis of 
SecHH #o(expt) is estimated separately. The estimation 
procedure is general and is outlined for Cu in Appendix 
iBl The agreement between theory and experiment in 
Table Hfl is certainly better than we should expect. 

The application of DFT to liquid dynamics research is 
currently making notable progress. From ab initio MD, 
melting properties of Si have been calculated by Sugino 
and Car (43|, and Wang et al. [3] calculated the carbon 
phase diagram. This work points to the possibility of an 
ab initio test of Eq. (fTUf for anomalous melting. High 
pressure melting curves have been calculated for Pb by 
Cricchio et al. [45| and also for Ta by Taioli et al. . In 
each case the shape of T m (P) suggests normal melting, 
and this again can be tested by calculating the quantities 
in Eq. (JTTJJ) . Kresse's summary of DFT calculations of the 
static structure factor and pair distribution function for 
group IIIB-VIB elements shows better results for metallic 
than nonmetallic liquids [47[ ■ While the metallic liquids 
should be well described by random valleys, the molec- 
ular character of As, Se, and Te, possessing strong and 
weak bonds, poses a challenge regarding the underlying 



5 



potential energy surface. A similar challenge is posed by 
Ge, whose ab initio static and dynamic structure factors 
compare well with experiment as shown by Chai et al. 
[48| , and whose primarily metallic liquid appears to have 
some tetrahedral coordination in its fluctuation spectrum 
at T m . 

The example of highly compressed Na has attracted 
much attention. It was predicted from theory by Neaton 
and Ashcroft [49[ that under compression crystalline Na 
will transform to low symmetry structures that include 
semimetallic behavior, and tend ultimately to semicon- 
ducting. Experiments by Hanfland et al. [501 ] and Syassen 
[5l| confirmed the structural changes to 120 GPa, and 
work by Gregoryanz et al. [52j revealed a change from 
normal to anomalous melting. This was in turn con- 
firmed by ab initio MD calculations by Raty et al. [53|, 
which also showed that the liquid undergoes electronic 
structure changes analogous to those in the solid. We 
notice that the anomalous melting regime in Na involves 
2s and 2p electrons entering the valence, and is related 
theoretically to the sequence of anomalous melting ele- 
ments Sn, Ga, Sb, and Bi mentioned in Sec. IIII Al 

A novel theoretical technique uses Monte Carlo pertur- 
bation theory to make accurate first principles calcula- 
tions of the liquid free energy at arbitrary temperatures 
[Hil I55I ]. This technique, as well as ab initio MD, will 
make possible more accurate tests of V-T theory than 
we can obtain through analysis of experimental data, as 
in the present study. 

IV. DISCUSSION AND CONCLUSIONS 

In this work, we are able to improve the formulation 
of the transit contribution to thermodynamic functions, 
and to carry out an ab initio test of the V-T theory of 
liquid entropy. The broad implications of these results 
will be discussed. 



The shape of the transit entropy curve in Fig. [T] can 
be understood from qualitative properties of transits. At 
very low temperatures (T « T m ) the system freezes into 
a single random valley and becomes an amorphous solid. 
The motion is entirely vibrational, there are no transits, 
and Str = 0, as is the corresponding Utr- Upon warm- 
ing, transits begin at a given temperature and both St r 
and Utr increase from zero. This behavior is seen in MD 
simulations of supercooled liquid Na, where the mean po- 
tential energy increases from its pure vibrational value of 
(3/2)Afc B T at around 140 K (Q, Fig. 4). This is con- 
firmed as a transit effect by observing that self diffusion 
increases from zero at approximately the same temper- 
ature ([13], Fig. 10). With increasing temperature, S tr 
and Utr saturate and have zero slope at the common 
temperature 6 tr - The reason for the saturation of Str 
and Utr, and their subsequent decrease with increasing 
T, is the truncation of the random valley potential sur- 
face at the intervalley intersections. A model for part 
of this decrease, the boundary effect, has been applied 
to the high-temperature specific heat of Hg 56]. Now 
we have more information than specific heat data, since 
Fig. Q] shows both the low-T increase and high-T decrease 
of the transit entropy. Figure [1] explains the previously 
disorganized behavior of liquid specific heat curves ( [56j , 
Fig. 2). Figure Q] will be helpful in modeling the statisti- 
cal mechanics of transits, and especially in modeling the 
transit free energy. 

In the past we have used T m as a scaling temperature 
for liquid properties. This is not satisfactory in principle, 
because T m depends partially on properties of the crys- 
tal. Not surprisingly, T m utterly fails as a liquid scal- 
ing temperature for anomalous melting elements ([l9j|. 
Figs. 3 and 4). We now have a two-component theory for 
the nuclear motion entropy, where each component has 
its independent scaling temperature, 8 Q (V) for S V ib and 
@tr(V) for Str- This is a purely liquid theory, with no 
parametric dependence on crystal properties. 



A. Transit contribution to thermodynamics 

In the original formulation, transits are accounted for 
through the multiplicative factor exp(A^lnui) in the par- 
tition function, where law = 0.8 is calibrated from en- 
tropy of melting data. Hence the transit free energy is 
—NksThiw, the transit entropy is the universal con- 
stant Nks In w, still a good approximation, but the tran- 
sit internal energy is zero. Now, in Sec. HI] the ex- 
perimental high-temperature entropy data are analyzed, 
and Fig. [T] is found for St r (V, T) at constant V. In 
principle one can integrate the constant-volume relation 
dU = TdS to find Ut r (V,T) up to a volume-dependent 
constant of integration. We then have the transit free 
energy F tr (V,T) to replace our original —NkBTlnw in 
the liquid free energy. This accomplishes the goal set out 
in the Introduction, and provides in principle the transit 
contribution to every thermodynamic function. 



B. Ab initio testing of V-T theory 

At the time of the original formulation no po- 
tential energy property of any random valley had been 
calculated. Indeed, the existence of the random and 
symmetric classes of valleys, and their contrasting po- 
tential energy properties, was only hypothesized in that 
first paper. At the time, we adopted the approxima- 
tion 6{){V m ) « 0Q{V m ) for normal melting elements, since 
9q was available from force-constant models calibrated 
to experimental dispersion curves [5 71 ] - This is still a re- 
spectable approximation, probably accurate to 3 — 4% on 
average for normal melting elements (but not for anoma- 
lous melting elements; see 58], Tables I and II). But 
now, with ab initio values of 9q(V), the correct theoret- 
ical S V ib(V,T) can be calculated from Eq. |2]) and com- 
pared with S V ii)(V, T) extracted from experiment. This is 
essentially the comparison made in Sec. IIIIBI But that 
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comparison goes much deeper than a casual glance would 
suggest. That test provides the following extensive sup- 
port of V-T theory. 

1. The test verifies the original hypothesis that the 
random valleys are numerically dominant and 
hence account for the entire statistical mechanics 
as N — > co, and that the random valleys are uni- 
form in their potential properties so that a single 
example is sufficient as N — > co. The verification 
results from the fact that a single random valley is 
used for the calculation of 0o for each liquid, while 
the experiment samples enormous numbers of val- 
leys of all types. 

2. The test verifies that vibrational motion in a har- 
monically extended random valley is the correct 
theory for that part of the experimental entropy 
which is identified with S v ib- This is because that 
part of the experimental entropy is reproduced by 
ab initio evaluation of Eq. ((2]) over the entire tem- 
perature range of the available experimental data 
for each liquid tested. 

3. The test is consistent with x(Vm) — 0.8-/Vfce for 
Na and Cu. The expected material dependence of 
x(^rn) i s a refinement remaining for future work. 

V-T theory is unique in that it offers a Hamiltonian the- 
ory capable of unifying equilibrium and nonequilibrium 
theories of liquid dynamics. Further testing will help 
to develop a robust theory. Techniques currently being 
applied in ab initio simulations can be employed to ad- 
vantage in the development of V-T theory 

APPENDIX A: ERROR ESTIMATES 

Error in S exp t{V m ,T) arises from experimental error 
in the high-T entropy data, and the experimental er- 
ror in the volume correction. We estimate the mean of 
each error over the range from T rn to T/j, and add their 
magnitudes to obtain the total mean experimental er- 
ror (5 S eX pt), as listed in Table HI Much of this error is 
already present in the data for S exp t(V m ,T m ). The re- 
maining error is T-dependent and has an average around 
zero for each liquid. This error is large enough to cause 
a noticeable error in the shape of Fig. [TJ 

The relative experimental error in 0o is <50o(expt)/0o, 
and is almost entirely due to (8 S expt)- From Eq. @ it 
follows that 

(56> (expt) _ {8 S expt ) . 
0o ~~ 3Nk B ' 



This gives a range of 1 — 5% for the experimental error 
in 0o in Table HI The error in 9 tr is entirely due to the 
T-dependent error, because 9t r depends on the shape of 
the curve and not on its magnitude. The error in 9 tr in 
Table U can reach 10%. 

Moments of the frequency distribution can be calcu- 
lated from DFT to an accuracy of 1% for elemental crys- 
tals. These crystal calculations are done in the infinite 
lattice model [2| , where an arbitrary number of Brillouin- 
zone fc-points is possible. In contrast, the liquid system 
has only 3iV normal models for an TV-atom system. Be- 
cause of the small system size (N — 150) for our DFT 
calculations, the present error in 9$(DFT) is allowed as 
3%. 



APPENDIX B: ESTIMATION OF SCALING 
TEMPERATURES FOR Cu 

The slope of St r {T m ) is Ct r (T m ), the transit contribu- 
tion to the specific heat, which can be found from data 
as follows. The experimental specific heat at constant 
volume is C eX pt, while the vibrational contribution in 
classical statistical mechanics is C v u, — 3Nk B , so that 
C tr — C exp t — 3Nk B — Cel. Matching this quantity, 
with error estimates, to the slope of a smooth curve 
fitted to Fig. [T] yields T m /9 tr in the range 0.8 - 1.3 
for Cu at T m . Then Fig. ED implies that S tr (V m ,T m ) 
is (0.78 — 0.80)Nk B . From experimental entropy data 
for Cu, S nuc (T m ) = 10.00Nk B , so that S mb (T m ) = 
(9.21 ± 0.01)Nk B , giving O = (171.4 ± 0.6) K. For Cu, 
(SSexpt) is just the estimated experimental entropy error 
at T m , namely 0.086 Nk B . 

This method becomes inaccurate when \Ct r \ is large. 
For a compilation of C v n, + Ctr for elemental liquids at 
melt, see [J]- 



ACKNOWLEDGMENTS 

For many helpful discussions and insights we thank 
Brad Clements, Giulia De Lorenzi-Venneri, Carl Greeff, 
Erik Holmstrom, and Travis Peery. This work was sup- 
ported by the U. S. DOE under Contract No. DE-AC52- 
06NA25396. 



[1] L. Boltzmann, Wien. Ber. 66, 275 (1872). Lattices (Clarendon Press, Oxford, 1954). 

[2] M. Born and K. Huang, Dynamical Theory of Crystal [3] M. Born and Th. von Karmen, Phys. Z. 13, 297 (1912). 



7 



D. C. Wallace, Phys. Rev. E 56, 4179 (1997). 
Macroscopic equivalence applies to one liquid at one den- [31 
sity. Statistical properties of random valleys vary with 
density, and vary from one liquid to another. [32 
D. C. Wallace, Statistical Physics of Crystals and Liquids 
(World Scientific, New Jersey, 2002). [33 
G. De Lorenzi-Venneri, E. D. Chisolm, and D. C. Wal- 
lace, Phys. Rev. E 78, 041205 (2008). [34 
W. A. Harrison, Pseudopotentials in the Theory of Metals 
(W. A. Benjamin, New York, 1966). [35 
N. W. Ashcroft, Phys. Lett. 23, 48 (1966). 
T. E. Faber, Theory of Liquid Metals (Cambridge Uni- 
versity Press, Cambridge, 1972). 

N. W. Ashcroft and D. Stroud, Solid State Phys. 33, 1 
(1978). [36 
N. H. March, Liquid Metals: Concepts and Theory (Cam- 
bridge University Press, Cambridge, England, 1990). [37 

D. C. Wallace and B. E. Clements, Phys. Rev. E 59, 2942 
(1999). [38 
B. E. Clements and D. C. Wallace, Phys. Rev. E 59, 2955 
(1999). [39 

E. D. Chisolm, B. E. Clements, and D. C. Wallace, Phys. 
Rev. E 63, 031204 (2001). [40 
G. De Lorenzi-Venneri and D. C. Wallace, Phys. Rev. E 

76, 041203 (2007). [41 
N. Bock, T. Peery, E. D. Chisolm, G. De Lorenzi-Venneri, 

D. C. Wallace, E. Holmstrom, and R. Lizarraga, Bull. [42 
Am. Phys. Soc. 53(2), J9:00004 (2008). [43 

E. Holmstrom, N. Bock, T. Peery, R. Lizarraga, G. De [44 
Lorenzi-Venneri, E. D. Chisolm, and D. C. Wallace, un- 
published. [45 

D. C. Wallace, Proc. R. Soc. London, Ser. A 433, 615 
(1991). [46 

E. I. Gol'tsova, High Temp. 4, 348 (1966). 
R. Hultgren, P. D. Desai, D. T. Hawkins, M. Gleiser, [47 
K. K. Kelley, and D. D. Wagman, Selected Values of the [48 
Thermodynamic Properties of the Elements (ASM, Met- 
als Park, OH, 1973). [49 
G. H. Shaw and D. A. Caldwell, Phys. Rev. B 32, 7937 
(1985). [50 
Yu. S. Trelin, I. N. Vasil'ev, V. B. Proskurin, and T. A. 
Tsyganova, High Temp. 4, 352 (1966). [51 
M. W. Chase, Jr., C. A. Davies, J. R. Downey, Jr., D. J. 
Frurip, R. A. McDonald, and A. N. Syverud, J. Phys. 
Chem. Ref. Data Suppl. No. 1 14, 61 (1985). [52 
E. E. Shpil'rain, K. A. Yakimovich, V. A. Fomin, S. N. 
Skovorodjko, and A. G. Mozgovoi, in Handbook of Ther- 
modynamic and Transport Properties of Alkali Metals, [53 
edited by R. W. Ohse (Blackwell, London, 1985), p. 435. 

N. B. Vargaftik, V. F. Kozhevnikov, V. G. Stepanov, [54 
V. A. Alekseev, and Y. F. Ryzhkov, in Seventh Sympo- 
sium on Thermophysical Properties, edited by A. Cezair- 
liyan (ASME, New York, 1977), p. 926. [55 
E. A. Brandes, Smithells Metals Reference Book (Butter- [56 
worths, London, 1983). [57 
R. R. Miller, in Liquid Metals Handboook, 2nd ed., edited 
by R. N. Lyon (U. S. Government Printing Office, Wash- 
ington, D. C, 1952), p. 38. [58] 
Y. Tsu, H. Suenaga, K. Takano, and Y. Shiraishi, Trans. 
Japn. Inst. Metals 23, 1 (1982). 

M. B. Gitis and I. G. Mikhailov, Sov. Phys. Acoust. 11, 



372 (1966). 

D. P. Almond and S. Blairs, J. Chem. Thermodyn. 12, 
1105 (1980). 

H. A. Spetzler, M. D. Myer, and T. Chan, High Temp. - 
High Press. 7, 481 (1975). 

H. Koster, F. Hensel, and E. U. Franck, Ber. Bunsenges. 
Phys. Chem. 74, 43 (1970). 

I. N. Makarenko, A. M. Nikolaenko, and S. M. Stishov, 
Inst. Phys. Conf. Ser. 30, 79 (1977). 

In the years between [Ts| and we realized that Tl 
should be removed from the analysis because available 
crystal data are for hep, while melting proceeds from 
bcc. This changes the average AS/Nks insignificantly, 
from 0.79 to 0.80. 

V. M. Glazov, S. N. Chizhevskaya, and N. N. Glagoleva, 
Liquid Semiconductors (Plenum, New York, 1969). 
D. A. Young, Phase Diagrams of the Elements (Univer- 
sity of California Press, Berkeley, CA, 1976). 
A. Jayaraman, W. Klement, Jr., and G. C. Kennedy, 
Phys. Rev. 130, 540 (1963). 

A. Jayaraman, R. C. Newton, and J. M. McDonough, 
Phys. Rev. 159, 527 (1967). 

J. Yamashita and S. Asano, J. Phys. Soc. Jpn. 29, 264 
(1970). 

D. B. McWhan, G. Parisot, and D. Bloch, J. Phys. F 4, 
L69 (1974). 

http://cms.mpi.univie.ac.at/vasp/ 

O. Sugino and R. Car, Phys. Rev. Lett. 74, 1823 (1995). 
X. Wang, S. Scandolo, and R. Car, Phys. Rev. Lett. 95, 
185701 (2005). 

F. Cricchio, A. B. Belonoshko, L. Burakovsky, D. L. Pre- 
ston, and R. Ahuja, Phys. Rev. B 73, 140103(R) (2006). 
S. Taioli, C. Cazorla, M. J. Gillan, and D. Alfe, Phys. 
Rev. B 75, 214103 (2007). 

G. Kresse, J. Non-Cryst. Solids 312-314, 52 (2002). 
J.-D. Chai, D. Stroud, J. Hafner, and G. Kresse, Phys. 
Rev. B 67, 104205 (2003). 

J. B. Neaton and N. W. Ashcroft, Phys. Rev. Lett. 86, 
2830 (2001). 

M. Hanfland, I. Loa, and K. Syassen, Phys. Rev. B 65, 
184109 (2002). 

K. Syassen, in High Pressure Phenomena, edited by R. J. 
Hemley, G. Chiarotti, M. Bernasconi, and L. Ulivi (IOS, 
Amsterdam, 2002), p. 251. 

E. Gregoryanz, O. Degtyareva, M. Somayazulu, R. J. 
Hemley, and H. K. Mao', Phys. Rev. Lett. 94, 185502 
(2005). 

J.-Y. Raty, E. Schwegler, and S. A. Bonev, Nature (Lon- 
don) 449, 448 (2007). 

C. W. Greeff and R. Lizarraga, Shock Compression of 
Condensed Matter - 2007, AIP Conf. Proc. No. 955 (AIP, 
New York, 2007), p. 43. 

C. W. Greeff, J. Chem. Phys. 128, 184104 (2008). 

D. C. Wallace, Phys. Rev. E 57, 1717 (1998). 

H. Schober and P. H. Dederichs, in Landoldt-Bornstein 
New Series Vol. 13a, edited by K.-H. Hellwege (Springer, 
Berlin, 1981). 

E. D. Chisolm and D. C. Wallace, Phys. Rev. E 69, 
031204 (2004). 



